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Abstract 

The ALICE Collaboration has measured the inclusive production of muons from heavy flavour de- 
cays at forward rapidity, 2.5 < y < 4, in pp and Pb-Pb collisions at -^/snn = 2.76 TeV. The p t - 
differential inclusive cross section of muons from heavy flavour decays in pp collisions is compared 
to perturbative QCD calculations. The nuclear modification factor is studied as a function of p t and 
collision centrality. A weak suppression is measured in peripheral collisions. In the most central 
collisions, a suppression of a factor of about 3-4 is observed in 6 < p t < 10 GeV/c. The suppression 
shows no significant p t dependence. 
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The study of ultra-relativistic heavy ion collisions is aimed at investigating the properties of strongly- 
interacting matter in the extreme conditions of high temperature and energy density expected to be 
reached. Quantum Chromodynamics (QCD) calculations on the lattice predict under such conditions 
the formation of a deconfined partonic phase, the Quark-Gluon Plasma (QGP), and chiral symmetry 
is restored [1 1. Heavy quarks (charm and beauty), abundantly produced at the Large Hadron Collider 
(LHC), are sensitive probes of the properties of the QGP Due to their large masses, they are created 
mainly in hard scattering processes during the early stage of the collision and subsequently interact with 
the hot and dense medium. In particular, measurement of open heavy flavour hadrons may probe the 
energy density of the system through the mechanism of in-medium energy loss of heavy quarks. The in- 
medium effects are usually quantified by means of the nuclear modification factor Raa of the transverse 
momentum (p t ) distribution. Using the nuclear overlap function from the Glauber model [2], Raa can 
be expressed as: 

1 dN AA /dp t 

R AA(jh) = -pfi—r--T—T7—, (1) 

where (7aa) is the average nuclear overlap function in a given centrality class. The term &N AA /&pi is the 
/^-differential yield in nucleus-nucleus (AA) collisions, while da pp /dp t is the /^-differential inclusive 
cross section in pp collisions. The value of Raa is unity for hard probes if no nuclear modification is 
present. A Raa value smaller than unity can arise from partonic energy loss as well as other nuclear 
effects. According to QCD, the radiative energy loss of gluons should be larger than that of quarks 
and, due to the dead cone effect Br|5)> heavy quark energy loss should be further reduced with respect 
to that of light quarks. The contribution from other interaction mechanisms, for instance collisional 
energy loss (7J[8l, in-medium fragmentation, recombination and coalescence l9TfTTl. could also lead to a 
modification of heavy-flavour hadron p t distributions in AA collisions. Finally, initial state effects lPT2l 
[T3l could complicate the interpretation of any deviation from unity of the Raa in terms of energy loss 
effects, particularly in the low p t region. The study of p-A collisions is required to quantify the role 
of initial state effects. The PHENIX and STAR Collaborations have reported a strong suppression of 
electrons from heavy flavour decays at mid-rapidity, in central Au-Au collisions at -^snn = 200 GeV at 
RHIC fl4l4"l6l . Recently, a significant suppression of D mesons ifTTl and J/i/^'s from B decays lfl8l was 
measured at mid-rapidity in central Pb-Pb collisions at ^/snn = 2-76 TeV by ALICE and CMS at the 
LHC, respectively. A complementary measurement of heavy flavour suppression at forward rapidity, at 
the same energy, is of great interest in order to provide new constraints on models which aim at describing 
the nuclear modification factor as partonic energy loss. 

In this Letter, we report the first measurement at the LHC of the production of muons from heavy flavour 
decays at forward rapidity (2.5 < y < 4), with the ALICE experiment [191, in PP and Pb-Pb collisions 
at y'j'NN = 2-76 TeV. The measured /^-differential inclusive cross section of muons from heavy flavour 
decays in pp collisions at yfs = 2.76 TeV is compared to perturbative QCD (pQCD) calculations. In- 
medium effects are investigated by means of the nuclear modification factor as a function of p t in 4 < 
Pt < 10 GeV/c, and as a function of collision centrality in 6 < p t < 10 GeV/c. 

The ALICE experiment is described in detail in [19|. The apparatus is composed of a central barrel 
(pseudo-rapidity coverage \r}\ < 0.9), a muon spectrometer (—4 < r\ < — 2.5 lJ) and a set of detectors 
for global collision characterization and triggering located in the forward and backward pseudo-rapidity 
regions. The VZERO, made of two scintillator arrays covering the 2.8 < T] < 5.1 and —3.7 < T] < 
— 1.7, is used for triggering, centrality determination and background removal. The two Zero Degree 
Calorimeters (ZDC), located at ± 114 m from the interaction point, are used in offline rejection of 
background events. The Silicon Pixel Detector (SPD), a two-layer central barrel that constitutes the 
innermost part of the Inner Tracking System (ITS), is included in the trigger logic. The SPD provides 
also the interaction vertex reconstruction. The muon spectrometer consists of a 10 interaction length 



'in the ALICE reference frame, the muon spectrometer covers a negative T] range and consequently a negative y range. The 
results are presented with a positive y notation. 
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(Ai) passive front absorber, a beam shield, an iron wall, a 3 T-m dipole magnet and a set of tracking 
and trigger chambers. Tracking is performed by means of five stations of cathode pad chambers, with 
the third station inside the dipole magnet. The tracking system is supplemented by two trigger stations 
of resistive plate chambers, behind a 1.2 m thick iron wall with thickness 7.2 X\. The latter absorbs 
hadrons that punch through the front absorber, as well as secondary hadrons produced inside it and low 
momentum muons, mainly from pion and kaon decays. 

The Pb-Pb data were collected during the 2010 run. The rate of hadronic collisions was about 100 Hz, 
corresponding to a luminosity of 1.3 • 10 25 cm _2 s _1 . The results presented in this Letter are based on 
the analysis of minimum bias (MB) trigger events. The MB trigger required the following conditions: 
a signal in at least two pixel chips in the outer layer of the SPD and a signal on each VZERO detector. 
The beam-induced background was reduced by using the timing information from the VZERO and ZDC 
detectors, and by exploiting the correlation between the number of hits and track segments in the SPD. 
Moreover, a minimal energy deposit in the ZDC was required in order to reject electromagnetic interac- 
tions. Finally, only events with an interaction vertex within ±10 cm from the center of the detector along 
the beam line were analyzed. Pb-Pb collisions were classified according to their degree of centrality by 
means of the sum of amplitudes of the signals in the VZERO detectors, as described in rt20ll2TI . The 
analysis was limited to the 80% most central events for which the MB trigger was fully efficient. This 
leads to a data sample of 16.6 • 10 6 Pb-Pb collisions which, in the following, will be divided into five 
centrality classes: 0-10%, 10-20%, 20-40%, 40-60% and 60-80% (the two last bins will be grouped 
together for the study of Raa(Pi))- The corresponding integrated luminosity is L- mt = 2.71 ±0.09 /lb -1 . 
The values of the mean number of participating nucleons and mean nuclear overlap function are given in 
TableQ] They were determined with the Glauber Monte-Carlo simulation assuming an inelastic nucleon- 
nucleon cross section of 64 mb [20]. The strategy of cuts applied to reconstructed tracks is similar to the 

Table 1: Mean number of participating nucleons ((N vart )) and mean nuclear overlap function ((7Xa)) for different 
centrality classes, expressed in percentiles of the hadronic Pb-Pb cross section. 



Centrality 


(AW) 


(Zaa) mb" 1 


0-10% 


357 ±4 


23.48 ± 0.97 


10-20% 


261 ±4 


14.43 ± 0.57 


20-40% 


157 ±3 


6.85 ± 0.28 


40-60% 


69 ±2 


2.00 ±0.11 


60-80% 


23± 1 


0.42 ± 0.03 


40-80% 


46 ± 2 


1.20 ± 0.07 



one used for pp collisions ll22l . Various selection cuts were used in order to improve the purity of the data 
sample. Tracks were required to be reconstructed in the geometrical acceptance of the muon spectrom- 
eter. A track candidate measured in the muon tracking chambers was then required to be matched with 
the corresponding track measured in the trigger chambers. This results in a very effective rejection of the 
hadronic background that is absorbed in the iron wall. Furthermore, the correlation between momentum 
and Distance of Closest Approach (distance between the extrapolated muon track and the interaction 
vertex, in the plane perpendicular to the beam direction and containing the vertex) was used to remove 
the remaining beam-induced background tracks that do not point to the interaction vertex, and fake tracks 
(tracks not associated to one single particle crossing the spectrometer). After these selections, the data 
sample consists of 10 • 10 6 muon candidates. The /?aa measurement of muons from heavy flavour de- 
cays will be performed at high p t (p t > 4-6 GeV/c) where the main background component consists of 
muons from primary pion and kaon decays. The Pb-Pb distributions are corrected for acceptance and 
for tracking and trigger efficiency (A x e) using the procedure described in ||22| . The global A x e is 
close to 80% for p t > 4 GeV/c. The dependence of the trigger and tracking efficiency on the detector 
occupancy, which is correlated with the collision centrality, was evaluated by means of the embedding 
procedure [23 1 . A decrease of the efficiency of about 4% ± 1% is observed in the 10% most central 
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collisions. 

The /?aa of muons from heavy flavour decays in the forward rapidity region is calculated according to 
Eq. (Q]), which can be written as: 



R 



-HF 



AA 



(A 



1 



d<Pb/d/V 



dAr p b pb 



/dp, 



(?aa) 



d< p 



-HF 



/dpt 



(2) 



where aWp bPb 



/dp t and dNp bpb ' /dp t are the inclusive muon and charged pion and kaon decay muon 
p t distributions at forward rapidity in Pb-Pb collisions, respectively. 



.j± , ]-[p 

The pp reference, dapp /d/? t , was obtained from the analysis of muon-triggered events collected 
during a pp run at y/s = 2.76 TeV, in March 2011, with integrated luminosity of 19 nb _1 after event 
selection cuts. The analysis technique from the event and track selection to the normalization, is the same 
as that described in ll22l . Figure Q] shows the measured /^-differential inclusive cross section of muons 



..a L ALICE pp ^s=2.76 TeV, li^HF in 2.5<y<4 
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\ 
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1 .9% normalization uncertainty not included 
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Fig. 1: (Color online) Transverse momentum differential inclusive cross section of muons from heavy flavour 
decays in 2.5 < y < 4, in pp collisions at ^/s = 2.76 TeV. The error bars (open boxes) are the statistical (system- 
atic) uncertainties. The solid curve and the band show FONLL [24 25 1 calculations and theoretical uncertainties, 
respectively. The FONLL calculations are also reported for muons from charm (long dashed curves) and beauty 
(dashed curves) decays, separately. The lower panel shows the ratio between data and FONLL calculations. 



from heavy flavour decays in the kinematic region 2.5 < y < 4 and 2 < p t < 10 GeV/c. In the range p t > 4 
GeV/c (pi > 6 GeV/c), regions of interest for the R^ A (p t ) measurement, the contribution of muons 
from primary light hadron decays (mainly primary pion and kaon decays) that was subtracted amounts to 
about 19% (12%) of the total yield. The error bars are statistical uncertainties. The open boxes represent 
the systematic uncertainties varying from 15% to 24%, depending on p t . This includes the contributions 
from background subtraction (ranging from a maximum of about 24% at p t = 2 GeV/c to 14% at p t = 10 
GeV/c), detector response (3%) and residual mis-alignment of tracking chambers (1% x p t , in GeV/c). 
The systematic uncertainty on the minimum bias pp cross section (1.9%), used in the normalization, is not 
shown. The data are compared to Fixed Order Next-to-Leading Log (FONLL) pQCD predictions B24II25I 
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(curve, with shaded band for the uncertainty). The ratio between data and FONLL calculations is also 
shown. The measured /^-differential inclusive cross section of muons from heavy flavour decays is well 
reproduced by the calculations within experimental and theoretical uncertainties, although at the upper 
limit of the predictions. A similar behaviour was also reported in pp collisions at yfs = 7 TeV ll22ll . 
The contributions of muons from charm and beauty decays from the FONLL calculations are displayed 
separately in Fig. Q] According to these predictions, the component of muons from beauty decays exceeds 
that of muons from charm decays for p t > 6 GeV/c. 

The p t distribution of muons from heavy flavour decays in Pb-Pb collisions at forward rapidity is ob- 
tained by subtracting the muon background component (mainly muons from primary pion and kaon 
decays) from the corrected inclusive muon /^-differential distribution. The presence of unknown nuclear 
effects, in particular medium-induced partem energy loss at forward rapidity, prevents subtraction of this 
contribution by means of Monte-Carlo simulations, as was done in pp collisions [22]. Hence, the contri- 
bution of muons from primary and K 1 * 1 decays at forward rapidity in Pb-Pb collisions was estimated 
by extrapolating to forward rapidity (2.5 < y < 4) the p t distributions of pions and kaons measured at 
central rapidity < 0.8) in pp and Pb-Pb collisions GSI and generating the corresponding p t distri- 
butions of decay muons with a simulation of the decay kinematics and of the front absorber. For the 
rapidity extrapolation, it was assumed that the suppression of pions and kaons is independent of rapidity 
up to y = 4. This assumption is motivated by the observation, by the ATLAS Collaboration, that the 
central-to-peripheral nuclear modification factor of charged hadrons does not show any T] -dependence 
up to T] = 2.5 within uncertainties B71 . The systematic uncertainty introduced by this assumption was 
conservatively estimated by varying R A ^ K (Pt) f rom (full suppression) up to two times its value. The 
entire background-estimation procedure is detailed in the following. 

The p t distribution of pions and kaons at forward rapidity in Pb-Pb collisions in a given centrality range 
is expressed as: 

^Pbpf/dA = (T AA )-da^ K± /d Pt • [< A K± ( A )],=o. (3) 

The mid-rapidity pion and kaon p t distributions measured in pp collisions were extrapolated to forward 
rapidity using [28]: 

d 2 A^ ±,K± /dAdy = [d 2 ^ ± ' K± /d/7 t dy] y= o-exp(g^) ) (4) 

with (7 y = 3.18. The latter is the average of the values obtained with the PYTHIA [29] and PHOJET El 
event generators. 

Then, the muon p t distributions in 2.5 < y < 4 in pp and Pb-Pb collisions were obtained by means of 
fast simulations using the resultant pion and kaon p t distributions as input. The effect of the front ab- 
sorber was taken into account by considering only pions and kaons that decay before reaching a distance 
corresponding to one interaction length in the absorber. 

The input charged pion p t distributions were measured up to p t = 20 GeV/c for all centrality classes 
used in the analysis. The kaon p t distributions were determined only at low p t . Therefore, the Kg p t 
distributions measured up to 16 GeV/c were used, considering that N(K + ) + N(K~) = 2 • N(K<!). A 
further extrapolation up to 40 GeV/c, by means of a power law fit, was needed. In addition, the Kg p t 
distributions were measured only for the 0-5% and 60-80% centrality classes. As a consequence, the p t 
distributions of muons from pion and kaon decays at forward rapidity were determined only in these two 
centrality classes. For the other centrality classes used in this analysis (Table 1), the <3NjLJ^ a K /dp t 

distributions were obtained by scaling the RjJ^ 7 * (pt) with the double ratio l?^ - K (px)/R A ^~ (p t ) 
which was found to be the same in the 0-5% and 60-80% centrality classes, within a maximum variation 
of 9% included in the systematic uncertainty. 

This procedure allowed us to estimate dAfp bpb /dp t and then to deduce the nuclear modification of 
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muons from heavy flavour decays at forward rapidity according to Eq. ([2]). The background contribution 
to the muon p t distribution increases with decreasing p t . Hence, in order to limit the systematic uncer- 
tainty on its subtraction, Raa was computed for p t > 4 GeV/c where this component is 7% (11%) of the 
total muon yield in central (peripheral) collisions. 

The systematic uncertainties on the Raa of muons from heavy flavour decays originate from the pp 
reference, the corresponding Pb-Pb yields and the average nuclear overlap function. The systematic 
uncertainty on the pp reference, previously discussed, is about 15-17% for p t > 4 GeV/c. The systematic 
uncertainty on the yields of muons from heavy flavour decays in Pb-Pb includes contributions from: 

- the inclusive muon yields in Pb-Pb collisions, about 6-10%, containing the systematic uncertainty 
on the detector response (3.5%), the residual mis-alignment (1% x p t , in GeV/c) and the centrality 
dependence of the efficiency determined with the embedding procedure (1%); 

- the yields of muons from primary pion and kaon decays in pp collisions at forward rapidity, about 
17%, due to the systematic uncertainty on the input mid-rapidity distributions, the extrapolation 
procedure (a y parameter) and the absorber effect (pion and kaon mean free path in the absorber); 

- the RfcA^ 71 (Pt)i about 14-17%, due to the systematic uncertainty on the input mid-rapidity pion 
p t distributions; 

- the «^"* ± ' K?t (Pt) /^aa - ^ Ot) double ratio, up to 9% at p t = 10 GeV/c; 

- the unknown suppression at forward rapidity for muons from primary pion and kaon decays. As 
mentioned, a conservative systematic uncertainty was considered by varying R AA ,K (pt) from to 
two times its value, with the additional condition that the upper limit does not exceed unity. 

Finally, the systematic uncertainty on the normalization includes the 1.9% uncertainty on the minimum 
bias cross section measurement in pp collisions and the uncertainty of 4.3% (centrality class 0-10%) to 
7.3% (centrality class 60-80%) on (r AA )- 

Figure [2] presents the Raa of muons from heavy flavour decays in 2.5 < y < 4, as a function of p t 
in central (0-10%, left) and peripheral (40-80%, right) collisions. The vertical bars are the statistical 
uncertainties. The /? t -dependent systematic uncertainties are displayed by the open boxes and include all 
the contributions previously discussed, except the normalization uncertainty that is displayed at Raa = 1 ■ 
A larger suppression is observed in central collisions than in peripheral collisions, with no significant p t 
dependence within uncertainties. 

The centrality dependence of the Raa of muons from heavy flavour decays was studied in the range 
6 < p t < 10 GeV/c where the contribution of muons from B decays becomes dominant in pp collisions 
according to the central value of the FONLL calculations (Fig. [TJ. The analysis was carried out in five 
centrality classes from 0-10% to 60-80% (Table [l). The resulting Raa is displayed as a function of 
(Npart) in Fig. [3] The contribution to the total systematic uncertainty which is fully correlated between 
centrality classes (filled boxes), including the pp reference and normalization, is displayed separately 
from the remaining uncorrected systematic uncertainty (open boxes). The Raa of muons from heavy 
flavour decays at forward rapidity exhibits a strong suppression with increasing centrality, reaching a 
factor of about 3^1 in the 10% most central collisions. 

The ALICE Collaboration has measured the production of prompt D mesons in 2 < p t < 16 GeV/c at 
mid-rapidity (\y\ < 0.5) [ 17 ] and the CMS Collaboration reported on that of non-prompt J/y from beauty 
decays, in 6.5 < p t < 30 GeV/c and \y\ < 2.4 [18]. The corresponding suppression of D mesons and J/y 
from beauty decays in those studies is similar to that reported here for muons from heavy flavour decays, 
although in a different p t and rapidity region. 
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Fig. 2: (Color online) /?aa of muons from heavy flavour decays in 2.5 < y < 4 as a function of p t , in the 0-10% 
(left) and 40-80% (right) centrality classes, in Pb-Pb collisions at -^/snn = 2.76 TeV. Open (full) boxes represent 
the systematic (normalization) uncertainty. Horizontal bars show the bin widths. 
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Fig. 3: (Color online) Raa of muons from heavy flavour decays as a function of the mean number of participanting 
nucleons, in 2.5 < y < 4 and 6 < p t < 10 GeV/c. The horizontal bars indicate the uncertainty on (N part ). 



In conclusion, we have reported on the first measurement of the production of high-/? t muons from heavy 
flavour decays at forward rapidity, in pp and Pb-Pb collisions at ^Ann = 2.76 TeV with the ALICE 
detector. FONLL pQCD calculations describe well the pp data within experimental and theoretical un- 
certainties, with the data being close to the upper limit of the model predictions. The Raa of high p t 
muons from heavy flavour decays indicate a clear suppression increasing towards the most central colli- 
sions. The measured suppression is almost independent of p u in the region 4 < p t < 10 GeV/c. These 
results provide clear evidence for large in-medium effects for heavy quarks in central Pb-Pb collisions at 
V'snn = 2.76 TeV The forthcoming p-Pb collisions will complement these measurements, by providing 
insight into the possible contribution of initial nuclear matter effects, although those are expected to be 
less important in the high p t region studied here. 
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